The Axiverse Extended: Vacuum Destabilisation, Early Dark Energy and 

Cosmological Collapse 
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A model is presented in tlie pliilosopliy of the "String Axiverse" of [3| that incorporates a coupUng 
of ultralight axions to their corresponding moduli through the mass term. The light fields roll in 
their potentials at late times and contribute to the dark sector energy densities in the cosmological 
expansion. The addition of a coupling and extra field greatly enrich the possible phenomenology of 
the axiverse. There are a number of interesting phases where the axion and modulus components 
behave as Dark Matter or Dark Energy and can have considerable and distinct effects on the 
expansion history of the universe by modifying the equation of state in the past or causing possible 
future collapse of the universe. In future such a coupling may help to alleviate fine tuning problems 
for cosmological axions. We motivate and present the model, and briefly explore its cosmological 
consequences numerically. 



I. INTRODUCTION 
A. The Mass Scale Of String Axions 

It was proposed in jl| that a generic prediction of string 
theory is the existence of many Hght axions. The La- 
grangian of such an axion can be characterised using two 
parameters: the symmetry breaking scale, fa (also re- 
ferred to as the axion decay constant), and the overall 
scale of the potential, A, appearing in the effective four 
dimensional Lagrangian: 
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where U (9) is some periodic potential. Bringing the ki- 
netic term into canonical form we define the field — faO, 
with Lagrangian: 



(2) 



where Vax{4') is again a periodic potential. Expanding 
the potential to quadratic order we find that the mass is 
given by: 
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Hence, any axion is equivalently parameterised by its 
mass and symmetry breaking scale. The authors of [l| 
then go on to assert that, because of the scaling of rua and 
fa with the parameters in a generic string theory com- 
pactification, and the complexity of such compactifica- 
tions, that one should expect fa to remain roughly fixed 
at some high scale fa ~ 10^^ GeV, whilst nia should dis- 
tribute roughly evenly on a logarithmic mass scale all the 
way down to the Hubble scale today of Hq ^ 10^^'^ eV. 

The argument can be summarised as follows. Axions 
arise from the existence of closed cycles in the com- 
pact space: one axion for each. The symmetry break- 
ing scale, /a, scales inversely with the action, S, due to 



non-perturbative physics on the corresponding cycle: 
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The action then typically scales with the area of the cor- 
responding cycle, so that across all axions in a given com- 
pactification volume fa will not vary over many orders of 
magnitude. The standard arguments then lead to stringy 
values of fa ~ 10^^ GeV 0- The crucial observation 
that enables the assertion about the distribution of ax- 
ion masses to be made is the exponential dependence of 
the scale A on 5: 
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where ^ is a mass scale: in the axiverse it is the geomet- 
ric mean of the supersymmetry (SUSY) breaking scale 
and the string/Planck scale. Therefore, as the areas of 
cycles scale over the volume of the compact space, and 
given that even the simplest compactifications contain 
hundreds of closed cycles of different orders, we can ex- 
pect axion masses to distribute roughly evenly on a log 
scale. 

It should finally be noted that in the axiverse scenario 
the lightest axions are not the standard QCD axion 
however they do owe their existence to it. String axions 
can be removed from the spectrum of light fields by tree 
level liftings towards the string/Planck scale, however if 
it is to solve the strong CP problem the QCD axion must 
escape such liftings. It would then be considered anoma- 
lous if the QCD axion were the only axion to escape such 
a lifting and remain light. 



B. String Axions As Dark Matter 

Axions and other ultralight scalar fields are a well mo- 
tivated candidate to make up a proportion of the dark 
matter in our universe 043' ^^'^ have also been moti- 
vated in many works to solve the problem of dark energy 
(see, for example, 0, 111)- Studies of the QCD axion Q 
(see for a review of the strong CP problem and the 
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axion as a possible solution) as a dark matter candidate 
[10I4T7I often depended heavily on the thermal proper- 
ties of the axion, particularly its temperature dependent 
mass. However, the thermal properties are specific to the 
non-perturbative instanton physics of QCD. In contrast, 
when we study axions arising in a generic context from 
string theory the non-perturbative physics of the ax- 
ion potential could arise from a variety of sources, and as 
such it is common to work with a simplified form of the 
potential with no temperature dependence. Therefore, 
cosmological axions in such a scenario have relatively 
simple dynamics: they acquire some initial conditions 
through spontaneous symmetry breaking of the under- 
lying global Peccei-Quinn (PQ) symmetry at the scale 
fa, and the field remains frozen at this value until the 
mass (provided by the potential due to non-perturbative 
physics, which arises at the scale ^. In the QCD case 
this happens near Aqcd) overcomes the Hubble friction, 
at which point the axion begins rolling towards and os- 
cillating about its minimum in the potential. Such fields, 
which acquire only small masses due to explicit breaking 
of a symmetry due to non-perturbative effects (instan- 
tons breaking the remaining shift symmetry of the PQ 
axion i n Q CD) are known as pseudo-Nambu-Goldstone 
bosons [i^ , and are the common motivation for the study 
of ultralight fields in cosmology. 

As a cosmological fluid the axions are purely gravi- 
tationally coupled: whilst frozen they behave as a cos- 
mological constant, and then make a short transition to 
oscillatory behaviour before behaving as a pressureless 
cold dark matter (CDM) component. The fraction of the 
total energy density in axions resulting from this non- 
thermal production depends upon the initial value the 
field acquires after symmetry breaking and its distance 
from the potential minimum (the misalignment angle), 
with a negligible thermal component from the (Planck 
and symmetry breaking scale-) suppressed couplings of 
axions to the standard model via higher dimensional op- 
erators. In such a scenario the initial misalignment angle 
therefore often requires tuning to produce a cosmology 
consistent with observations if we assume the existence 
of an axion with a certain mass. The tuning may be to 
small values of the misalignment angle if we wish not to 
overdose the universe with heavier axions, or to large val- 
ues of the misalignment angle if we would like the dark 
matter to contain some significant fraction of ultralight 
axions @. It should be noted, however, that the measure 
of fine tuning in such scenarios depends not only on the 
potential, which is unknown and can in some cases be 
non-periodic (in which case the concept of tuning is 
more ill-defined), but also on the details of infiationary 
physics [^-i3. 

These fine tunings are a problem and a blessing for 
ultralight fields. In the "anthropic boundary" for ul- 
tralight axions was discussed. This refers to the region 
of axion parameter space where it is no longer possible 
for the axions to make up an 0{1) fraction of the dark 
matter, due to the periodic nature of the field making 



ip = fair the maximum possible misalignment. Unless the 
periodicity is broken as was done in Q and assumed in 
Q, or anharmonic effects in the potential are accounted 
for ^145, ultralight axions cannot make up large enough 
fractions of the dark matter to be readily observable, but 
consequently pose no danger of "overdosing" the universe 
and producing a Hubble rate inconsistent with observa- 
tions. Part of the motivation for this work is to provide 
a mechanism that decouples the axion density somewhat 
from its initial misalignment angle via a tracking mecha- 
nism [U, [25| , which may allow for observable densities in 
light axions, or to accommodate larger fa for the heavier 
axions. 

Recently a specific construction in string theory real- 
ising the "axiverse" of [l| was given [2^, [23] ■ The au- 
thors of [2^ [131 suggested that non-thermal processes 
are particularly important for the cosmology of such a 
model. In this paper I will propose a simple extension of 
the axiverse that includes further non-thermal processes 
that couple the axion fluid and greatly enrich the phe- 
nomenology of axion dark matter in a way analogous to 
the enrichment of dark energy phenomenolo gy a rrived at 
through the study of coupled quintessence [28|. In this 
framework the behaviour of the axion dark matter com- 
ponent is no longer a simple transition, and one hopes to 
look for novel features in the late time expansion history 
of the universe, as has been done fruitfully in many works 
on dark energy (see, for example, [29l. [SOj). 

The perturbations of ultralight dark matter axions and 
their effect on structure formation has been studied in 
many previous works [HI, [^ , the basic result being that 
the large Compton wavelength gives rise to a quantum 
pressure and suppresses structure formation below that 
scale 3 , in a way largely analogous to the presence of a 
massive neutrino [^ . The study of perturbations in this 
extended scenario will be left to a future work. 



C. Moduli 

Moduli are scalar fields in string theory that control 
the size and shape of the compact manifold. That there 
are many of these and they can contribute significantly 
to the cosmological energy density if they are allowed 
to roll in their potentials is known as the "cosmological 
moduli problem" , and moduli stabilisation is an impor- 
tant problem in string theory, related to the landscape, 
the existence of stable vacua, and the cosmological con- 
stant problem (see [sj and references therein). A typical 
potential for moduli stabilisation in the LARGE volume 
scenario for a modulus x is of the parametric form: 

where B, D and C are for our purposes free parameters, 
with the only restriction that for a stabilising minimum 
we have 2B > D. The typical scales of these parameters 
will be discussed shortly. We will also add an arbitrary 
(cosmological) constant, A (not to be confused with the 
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scale of the axion potential), to the potential Vmod to ac- 
count for the frozen/ vacuum expectation values of other 
moduli and standard model fields, the vacuum energy it- 
self, and the bare gravitational A in Einstein's equations, 
thus allowing for a late time de-Sitter expansion and ac- 
ceptable phenomenology. This is a prescription necessary 
in most all moduli stabilisation, since the vacuum found 
is always anti-de-sitter. 

The modulus dynamics are therefore expected to be 
a small perturbation about the standard model for dark 
energy, similar to models of varying equation of state, 
quintessence or Early Dark Energy (EDE) (For a com- 
prehensive review of dark energy theory and phenomenol- 
ogy, see d^l-). This prescription, like many cosmological 
models, is not answering the cosmological constant prob- 
lem, only giving it some alternative phenomenology. The 
fate of the universe given this scenario is discussed in Sec- 
tion |llll 

In the LARGE volume moduli stabilisation scenario 
of [m, [11] the lightest modulus is the volume modulus, 
which like the dilaton controls the overall scale of cou- 
plings in the model, but we can expect the presence of 
many light moduli. There are moduli corresponding to 
each closed cycle in the compact space, and as such there 
is a modulus for each axion. In particular there will be a 
modulus controlling the area of a cycle and as such: 

S-X (7) 

The modulus mass scale that leads to rolling in the poten- 
tial is often tied to the axion mass, so that the presence 
of cosmologically rolling axions suggests the presence of 
cosmologically rolling moduli. In this way we extend the 
axi verse. 

The displacement of moduli from their stabilised val- 
ues (vacuum destabilisation) by astrophysical processes is 
not a new idea: it has been exploited before in chameleon 
models (where destabilisation is attained via coupling 
to density) and other similar scenarios (see, for exam- 
ple |34| - [37j ). The variant on the scenario that is pro- 
posed here is as follows. The lightest axions and moduli 
have masses below the Hubble scale and are stabilised by 
Hubble friction, contributing to an effective cosmological 
constant. The heaviest moduli are stabilised by poten- 
tials of the form Eqn. [SJ which includes the moduli for 
standard model couplings ^ . There are then cosmological 
axions that roll in their potentials at late time and con- 
tribute to the dark matter density and perturbations in 
a distinct way from standard CDM. As such, some mod- 
uli may also roll in their potentials on cosmological time 
scales, which will lead to a rolling of the scale of the axion 
potential and consequently a rolling of the axion mass. 
It is this feature of coupling that we hope to exploit in 



looking for novel features in the cosmological expansion 
rate and equation of state in the dark sector. It is worth 
stressing again that cosmologically rolling moduli are not 
general, but we explore the phenomenology of allowing 
such a scenario for some of the moduli. 

A final word of warning should be made about the ad- 
dition of cosmologically relevant scalars to any model. 
Without a symmetry, such as the axion shift symmetry, 
which forbids scalar couplings directly to such terms as 
F^^F^'^ (although they might be induced at loop order 
like the axion-photon coupling), scalars will always ap- 
pear multiplying such terms and if they are light will 
induce long range, gravitational strength "fifth forces". 
This problem is generic to scalar quintessence models, 
though often ignored. The chameleon models are geared 
towards solving this problem [3^ . This work will not ad- 
dress such fifth force constraints, but they should be born 
in mind when making any detailed analysis, and will also 
be the subject of future investigation. 



D. Comments 

In this model we are working in the Leib- 
nizian/Panglossian philosophy of the authors of T] to use 
the vastness of the string landscape to look for general 
properties exploitable for cosmological phenomenology 
i.e. realistic string compactifications and scenarios for the 
moduli are complicated and varied but have many general 
and model independent features j33|. In the literature 
there has been a long standing link between string the- 
ory and inflation physics (see (39jand references therein) . 
In particular it was the aim of [39| to connect inflation- 
ary observables to topological properties of the compact 
space; this is analogous to "cohomologies from cosmol- 
ogy" in [l[ . However in the simplest axiverse scenario for 
dark matter axions there are two parameters per axion 
that can be constrained using cosmology: the axion mass, 
ma, and the axion fraction, / = f2a,/r2,„. The fraction / 
is determined by the initial misalignment angle and tells 
us nothing about the compact space or vacuum, so any 
cosmological bounds based on the expansion rate or for- 
mation of structure can only limit the contribution from 
a given axion or number of axions of given masses rather 
than place definite constraints on the number and size 
of closed cycles that determine the masses. There could 
be many light axions of many masses but if their con- 
tribution to the energy density is too small then we will 
not observe them ^. An observation of the effects of a 
number of axions, modulo a prior on the potential and 
degree of fine tuning, may also hope to place bounds on 



^ We will not investigate any evolution of the couplings with cosmic 
scale. The possibility and implication of such effects in this model 
will be the subject of future work 



^ The existence of an axion field coupling to standard model fields 
can be constrained in many other ways: for example, by light 
shining through a wall experiments (for a recent review , see 1401 ). 
and astrophysical processes such as those explored in |4ll444ll 
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their common axion decay constant, fa- Assumptions 
about inflation can also bound fa, as discussed in ^19r - 
[23| . and vice versa an assumed ultralight axion bounds 
Hi from considerations of isocurvature perturbations, as 
discussed in [l7| . 

Considering axions and moduli in the landscape and 
whether they survive to be observed is closely linked to 
the question of SUSY breaking, and is discussed in [45| . 
We will discuss it no further here, except to stress again 
that issues of naturalness and fine tuning within the land- 
scape will not concern us here: we simply motivate a the- 
oretically plausible, phenomenolgically viable, and obser- 
vationally testable scenario. 

Extending the axiverse scenario as proposed here gives 
much greater scope for direct connections between late 
time cosmology and beyond the Standard Model/string 
physics, akin to those already fruitfully explored in infla- 
tionary physics, but distinct from those connections al- 
ready explored in tackling the dark matter problem with 
weakly interacting massive particles (WIMPs) in the con- 
text of supersymmetric extensions of the standard model 
(the MSSM and its progeny) , or in directly addressing the 
cosmological constant problem. However, as phenome- 
nologists we must be careful about any statements we 
make about fundamental physics based on any results ob- 
tained using our models. In a parameterised phenomeno- 
logical model such as the one proposed here, the param- 
eters should be taken as just that: it is only in the con- 
text of a full string model such as those in f^, [2^ Isi], 
(with all the assumptions that go into constructing such 
a model) that the parameters take on their physical high 
energy physics meanings. 



II. THE MODEL: COUPLED AXIONS IN THE 
DARK SECTOR 

A. Equations of Motion 

The coupled axion-modulus Lagrangian takes the fol- 
lowing form: 

(8) 

where (f> is the axion, x is the modulus, B, D are dimen- 
sion four parameters for the modulus potential, C is re- 
lated to the overall volume of the compact space in string 
units, C is related to the instanton action, = fJ^ / fa 
and we make the simplifying assumption to work with 
only the mass term for the axion. In the interests of 
economy of parameters we will often take C = C. Also, 
when C ^ C the problem of minimising the potential 
becomes much less tractable generically, and for specific 
values of C, C many minima appear (e.g. 19 minima with 



(7 = 10,(7 = 1), which is related to the emergence of the 
landscape in string theory and the analysis of which is 
beyond the scope of this work. 

In this work we will only be concerned with the ho- 
mogeneous background fields, so will use the notation 
(j) = (f)Q = (j)o{T), where r is conformal time, and simi- 
larly for the modulus. 

The energy momentum tensor and equations of motion 
for the coupled system follow in the usual way from the 
Lagrangian. For a homogeneous Friedmann-Robertson- 
Walker metric, with scale factor a, in conformal time: 

fa 

X + 2HX- Ca^2Be-^''^ - De-^^) = ic^ft'e^c^^^^ 

2 fi 

(9) 

where over dots denote derivatives with respect to con- 
formal time, and T-L = a/ a. The energy momentum ten- 
sor for the combined axion-modulus system has the form 
of a perfect fluid with energy density p and pressure P: 
= -p, = P5' j. This gives: 

p=^(0' + x') + mx) 

(10) 

Due to the coupling, only these combined quantities obey 
the conservation equation p — —3'H{p + P). 

The form of the potential suggests a natural splitting 
of this into components due to the axion, subscript 0, 
and modulus, subscript x'- 
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(11) 

though in certain cases this distinction should be looked 
at more carefully [46| . 

With these definitions we will investigate the scalings 
of the energy density for axion and modulus compo- 
nents, and also the combined system. It is not only 
the scaling of the energy density that effects the cos- 
mological expansion history: we will also investigate the 
equation of state, given by: Wi = Pi/ Pi for components 
i = ax, mod, ax+mod, tot, which will effect the expansion 
rate in the usual way [U 13 ■ 
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The Friedmann equation is: 



parameters'^: 



SttG 



-a p 



(12) 



In addition to the axion-modulus system of Eqn. [10] we 
win consider components of the energy density coming 
from radiation, p^, CDM, pc, and a cosmological con- 
stant, Pa, all of which will redshift in the usual way. This 
now completes the description of the system. 



The Scales of Parameters 



This model is to be understood phenomenologically, 
and thus all the parameters will be taken as free when 
searching for interesting cosmological features, however 
it will be useful to have some idea of the natural scales 
in relation to the units used in numerical solution of the 
equations. Firstly, we scale the reduced Planck mass, 

SttG = TJ3-, out of the Friedmann equation by rescaling 

pi 

the fields to be in Planck units: (j)/Mpi, x xl^pi, 
and absorbing factors of M^; into B, D and the densities 
of the standard ACDM components. Next we change 
time variables to work in units of H^: t — >■ Hqt. This can 
be divided through the equations of motion and absorbed 
into the parameters and densities, so that the densities 
are now: p — > p/(HlM'^^), and Pi(a)/3 — ^i{a). Finally 
we express all the parameters in the potential in Planck 
units, natural for a string inspired model, as X ^ M^^X , 
with X the parameter and x its mass dimension. Thus, 
finally we have: 



B 



D 



M = 



Pa 



B 



D 



HI 



SUSY 



B < 10^^° 
D < IQi^o 
M < 10^^ 

For the fields, a large value of > 1 will be trans- 
planckian, and any value (f) > T^{fa/Mpi) will represent a 
departure from the periodic nature of the axion, whereas 
X is an area/volume, so a large value in Planck units is 
not problematic, and is in fact what one expects in a 
LARGE volume scenario for moduli stabilisation. The 
initial conditions on the fields are free parameters in the 
model: x representing a position in the landscape at the 
end of infiation, and (j) being selected by spontaneous 
symmetry breaking. The flatness condition fixes the den- 
sity in the axion-modulus system. Choosing the densities 
in A and standard CDM to be close to their observed 
values allows the axion-modulus system to be set as sub- 
dominant. Thus for now we ignore the question of fixing 
the initial conditions so as to obtain Hq at its observed 
value, since the subdominant components will not cause 
much variation away from this {Hq = 1 in our units). 
The appropriate initial conditions for a universe contain- 
ing a significant fraction in an axion-modulus component 
will be the subject of future work. 

We expect the dimensionless parameters C and C to 
both be 0(1). 



C. Comments 

At this stage some comments on the system of equa- 
tions in relation to other models in the literature will 
be useful. A brief comparison will be made to three 
models: [1, H^, , stating the main similarities and dif- 
ferences. The take-home message, though, will be that 
these models, whilst interesting, are connected to very 
different sources of new physics. They are for the most 
part motivated by scalar-tensor theories modifying the 
gravity sector to build dark energy models, whilst the 
model presented here is firmly cast in the context of HEP, 
strings, and the landscape, with dark matter candidates 
and modifications to the dark energy sector a by-product 
useful for phenomenology. 

In [4^ the scalar field analogous to our modulus is cou- 
pled non-minimally to gravity via the term ip^R, and has 
its own scalar potential. The field is also used to Higgs 
the Dirac dark matter particle via a Yukawa coupling 
tp f {(f>)'ip , thus introducing a cf) dependent mass to the 



All the parameters on the left hand side are now dimen- 
sionless, and it is these that will be used when quoting 
results. 



Using M2 /i/, 



10120^ /a ^ IQi^GcV, and Planck 
scale SUSY to give an upper bound on M = 1 gives 
an idea for the approximately natural scales of all the 



^ TeV scale SUSY is obviously more attractive from the point of 
view of resolving the hierarchy problem. However, B and D 
should also be set at the same magnitude, e.g. from potentials 
due to gaugino condensation, and our conclusions only rely on 
{/aVM2)M2~B~i3. 
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dark matter. However no explanation for this coupling is 
given in terms of fundamental particle physics (general 
couplings of this form can, however, be constructed in 
scalar-tensor theories by transforming between the Jor- 
dan and Einstein frames: see for example (sojV Assum- 
ing they have a standard neutralino dark matter particle, 
then the scalar field will be a Higgs of some supersymmet- 
ric extension to the Standard Model and one may worry 
about effects both on the freeze out mechanism for the 
neutralino, or on other effects where the true Higgs rolls 
and causes mass changes in the Standard Model parti- 
cles. This model building issue aside, the authors go on 
to explore the effects of different potentials and coupling 
terms that allow matching to the standard cosmology. 

There are some important distinctions between the 
model of [i^ and the one presented here, the first being 
that the axion-modulus model makes no change in the 
gravity sector, and as such the FRW equations take their 
standard form, though [i^ can be brought into stan- 
dard form via a conformal transformation: this is the 
usual degeneracy of Modified Gravity theories to scalar 
field theories of dark energy (again, see [50]). Secondly, 
the dark matter sector in |49j is, apart from the cou- 
pling, assumed standard, whereas the axiverse model is 
motivated by non-standard dark matter components. Fi- 
nally, many scalar potentials and couplings are explored 
in [49| within a general framework, whereas in the axion- 
modulus model the forms are essentially fixed via the 
motivation in string theory, and would be formally fixed 
in any specific string realisation of the model. The mo- 
tivation of [4^ is to explore late time dynamics of dark 
energy; the motivation here is to investigate dynamics in 
the matter sector and finds dark energy dynamics as a by 
product. Thus both models explore the same idea of dark 
sector couplings and non-standard dynamics, the model 
of 49] being more phenomenologically general for cosmol- 
ogists and dark energy phenomenologists, the model here 
being more theoretically and phenomenologically general 
from a HEP point of view. 

Coupling axions and moduli bears much similarity 
to "axion-dilaton cosmology" [1, [5l|. In these mod- 
els the axion has no potential of its own but is cou- 
pled in its kinetic term to the dilaton after a confor- 
mal transformation renders the gravity sector standard: 
L C —\e^''{dx)^, where a is the dilaton, and x the 
axion. The dilaton potential is then also exponential 
Ae""*"^, appearing on the cosmological constant term af- 
ter the conformal transformation. This model is then 
entirely constrained with only two free parameters. The 
dynamics in the scalar field sector consists of the exis- 
tence of attractors bringing the equation of state period- 
ically into accelerating phases, with possibly observable 
consequences as explored in [29]. The highly constrained 
axion-dilaton model should serve as a useful heuristic 
guide when thinking about the coupled scalar field dy- 
namics of the axion-modulus model, the important dif- 
ferences being that axion-modulus model has an axion 
mass term with a coupling on it, but standard kinetic 



term, and the potential for the modulus has a finite field, 
negative potential minimum, in contrast to the dilaton. 
We will expect this to make quantitative and qualitative 
changes to the equation of state evolution. 

Finally the model of 46:] bears the most resemblance 
to the axion-modulus model with the mapping of axion 
to "geon" X, and modulus x to "cosmon" kp. the only im- 
portant difference being that the modulus potential has 
a minimum, whilst the cosmon potential, being a dila- 
ton, does not. The analysis of [4g] will accordingly be 
very useful for helping us understand the axion-modulus 
system. However, I choose to include a standard cos- 
mological constant as well as the modulus, for the place 
holding reasons mentioned in Section HI and consider this 
good practice. Whether or not A can be removed from 
the axion-modulus system, with all the dark energy given 
from the scalar fields, will be the subject of future work, 
and if possible may have important consequences in rela- 
tion to string theory. I also choose to include a standard 
dark matter component because we know that if ultra- 
light scalar fields rolling on the time scales of interest 
here and in (46j suppress structure formation according 
to their fraction [^,0], then, just like massive neutrinos, 
they cannot make up all of the dark matter, as is assumed 
in [46]. 

This concludes the discussion of the model. 



III. RESULTS 

A. Example Cosmology 1 

This first example is concerned with parameter values 
close to those considered natural in Section|lIl The evolu- 
tion of the various components of the density p obtained 
from numerically integrating the equations of motion are 
plotted as a function of scale factor a in Figure[Tl with the 
associated fi's plotted in Figured] for parameter values: 
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= 10^2 
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= 10^20 
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= 10"9 
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= (7 = 10 
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= 1 


Xt 


= 25 




= 0.7 




= 0.2 




= 8 X 10" 



(13) 

In Figure [T] there are a number of qualitative features 
worthy of comment. Firstly, the logarithmic scale pre- 
vents us from showing the small negative energy density 
associated to the modulus at early and late times, we 
only see it emerge onto the plot between scale factors 
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FIG. 1: Densities as a function of scale factor for the param- 
eters of Equation 1131 Notice the tracking dynamics of the 
axion-modulus system between 10~^ % o, 1^ 10~^, when the 
modulus gains positive energy density, and the end of this 
tracking at a > 10~^ caused when the axion field begins os- 
cillations and the axion density makes its standard transition 
to CDM-like behaviour. 

10^® a < 10~^ as it enters an attractor scaling so- 
lution. During this time the modulus is displaced from 
its initial value and its equation of state becomes kinetic 
dominated, w — \^ the additional energy density being 
kinetic, as demonstrated in Figure [31 This transit causes 
a corresponding evolution of the axion mass, as demon- 
strated in Figure [HI It also induces a tracking behaviour 
in the axion density whilst the modulus is rolling. The 
axion field then begins its usual oscillatory behaviour and 
the axion behaves as dark matter for a > 10~^. 

The effects of this are best viewed in terms of the 
equations of state for the combined systems. The axion 
equation of state, Wax, and the axion-modulus equation 
of state, Wax+mod, are shown in Figure ID The axion 
equation of state differs from the usual case of a quick 
transition between w = ~1, and oscillations averaging to 
w — 0. It is in the combined equation of state that we 
see tracking behaviour as w tries to follow the equation 
of state of the dominant component, before the axion os- 
cillations begin, which are the cosmic trigger event that 
destroys tracking and restabilises the modulus (in this 
context, stabilisation is defined by Wmod ^1)- Once 
the axion field begins oscillations it causes the equation 
of state to oscillate and the axion pressure averages to 
zero, i.e. the axion behaves as pressureless dark mat- 
ter. The final cosmology today at a = 1 has a negligible 
negative component of energy density in the frozen mod- 
ulus, whilst an ultralight dark matter axion makes up a 
fraction of the critical density of order that in radiation. 

The resulting behaviour of the total equation of state 
for all fluids is plotted in Figure[Sl where we see that it de- 
velops a kink around scale factor a ~ 10~^ away from its 
usual ACDM evolution through matter-radiation equal- 
ity, caused by the presence of a significant component of 
the fractional density due to the axion-modulus system at 
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FIG. 2: Fraction of the critical density, f^, as a function of 
scale factor for the parameters of Equation 1131 
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FIG. 3: Evolution of the modulus equation of state ui as a 
function of scale factor for the parameters of Equation 1131 
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FIG. 4: Evolution of the axion and combined axion-modulus 
equation of state w as a function of scale factor for the pa- 
rameters of Eauation ll3l Notice the tracking behaviour of the 
combined equation of state, and the end of this when axion 
field oscillations cause oscillations in the pressure, averaging 
to zero. 
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FIG. 5: Evolution of the total equation of state w as a func- 
tion of scale factor for the parameters of Equation 1131 Note 
the appearance of a feature around a ~ 10"'' that departs 
from the standard smooth ACDM evolution. Compare this 
to Figure [2] where an overshoot as the axion- modulus system 
enters its tracking solution causes a significant contribution 
to the critical density at this scale factor. 




FIG. 6: Evolution of the axion mass as a function of scale 
factor for the parameters of Equation 1131 Compare this to 
Figure [S] the rolling occurs whilst Wmod ~ 1 and the modulus 
has kinetic energy. 



this time, as demonstrated in Figure [5J This is EDE-like 
behaviour, although the equation of state is not dragged 
low enough to cause an early period of acceleration. 

Having identified the main features of cosmology in 
this example, we now turn to briefly assess the depen- 
dency of these features on the parameters. A full descrip- 
tion of the system in this way with its various degenera- 
cies will require analysing it as a system of autonomous 
equations in the phase plane and the identification of the 
fixed points (25l [28j , which is left for a future work. 

Decreasing the axion initial value, 0;, can cause a sig- 
nificant change in the modulus behaviour. With C = C 
the condition for the modulus to have a finite real mini- 



mum is given by: 

2D > M^(j)^ 

For the values used in Equation [13] this can only oc- 
cur as the axion decays, and the = minimum is at 
X = 0.3. Now, the dependence on axion initial condi- 
tion can be most easily seen by plotting the potential 
V{(l),x)- In Figure [7] the potential is plotted for cj) in the 
range { — 1,1}, whilst in Figure [8] it is plotted with (j) in 
the range {—0.01, 0.01}. We see that, for small amplitude 
axion oscillations the exponential descent into the mod- 
ulus minimum for (j) — shows up strongly, whereas this 
feature is hidden on the same scale when the oscillations 
of the axion have a larger amplitude. We can similarly 
remove the axion from the spectrum almost entirely by 
increasing C by an order of magnitude, making the axion 
effectively massless and allowing the modulus to rapidly 
reach its minimum. 

Thus it is the axion oscillations which here stabilise the 
modulus away from its true minimum (as is visible com- 
paring Figures [T] and [3]) . If the modulus falls into this 
minimum then a large negative potential is generated 
causing the universe to collapse [52|- That this occurs 
prior to a = 1 for a small axion initial field value corre- 
sponding to reasonable initial misalignment 9i = 1 (not 
shown) rules out these particular parameters, and will 
require tuning in these models to avoid it. This decay 
and collapse may occur in the future as the amplitude 
of axion oscillations decays ^] for any model that looks 
viable today, and the parameters can then be used to 
estimate the lifetime of the universe [s^ . The analysis 
of a collapsing universe in this model will be the subject 
of future work. Some further comments on entry into a 
collapsing phase will be made in Section IIII Bl 

Even a sight decrease in C* to 8 or 9 also has a dramatic 
effect on the potential, allowing for the appearance of 
new minima as the axion field oscillates about zero, and 
the and D terms in the potential play off against 
one another. These new minima are sharp and highly 
localised in field space. Entry into them occurs as the 
modulus grows through its tracking solution and axion 
oscillations decay. A brief dip through them leads to a 
short period of negative potential domination, which the 
Friedmann equation is incapable of dealing with, and as 
such the analysis of this region of parameter space is left 
for a future work. 



B. Example Cosmology 2 

The purpose of this example is to demonstrate the 
large freedom in choosing values for parameters in these 
models by producing a similar and viable cosmology to 
that presented in Section IIII Al but with parameters 
many orders of magnitude different. This also shows that 
we expect many degeneracies in the parameters, with 
only the ratios of some being relevant. Specifically, the 



9 




of Equation 1131 

FIG. 9: Densities as a function of scale factor for the param- 
eters of Equation 1141 




FIG. 8: The potential V{(t), x) of Equation E] plotted for x 
in the range {24, 25}, </> in the range {-0.01, 0.01}, for the 
parameters of Equation 1 131 

parameters in this example are: 



M = 


10^ 




B = 


10*^ 




D = 


10^ 




C = 


c = 


= 10 




10^ 




Xi = 


10" 


5 


r!A = 


0.7 






0.2 






8 X 


10" 



(14) 

The evolution of the densities, fi's, equations of state, 
axion mass, and the axion field, are shown in Figures |9l 
[inillllllllllElandllSl 

This example shows best the possibility of the axion- 
modulus system to display tracking EDE behaviour. The 
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FIG. 10: Fraction of the critical density, fi, as a function of 
scale factor for the parameters of Equation 1141 
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FIG. 11: Evolution of the modulus equation of state to as a 
function of scale factor for the parameters of Equation 1141 



10 




1000.0002 



FIG. 12: Evolution of the axion and combined axion-modulus 
equation of state w as a function of scale factor for the pa- 
rameters of Equation 1141 
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FIG. 13: Evolution of the total equation of state w as a func- 
tion of scale factor for the parameters of Equation [141 
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FIG. 14: Evolution of the axion mass as a function of scale 
factor for the parameters of Equation 1141 
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FIG. 15: Evolution of the axion field as a function of scale 
factor for the parameters of Equation 1 141 



densities of both fields are always positive and show a 
scaling with the dominant component of energy density. 
Whilst the modulus equation of state is kinetic domi- 
nated once rolling begins, w = 1, the axion equation of 
state remains always potential dominated, w = —\, and 
the axion slowly rolls, never beginning oscillations as the 
mass is exponentially damped below the Hubble scale 
(see Figures [TH [T5|) . Tracking for the combined system 
thus persists into the present epoch, and should persist 
indefinitely, being absent the trigger event of axion oscil- 
lations to end it. In this case, the splitting of the energy 
density between the two components as done in Equa- 
tion [11] is not really so clear and it makes sense to speak 
more in terms if the combined axion-modulus system as 
a quintessence fluid. 

In this example the equation of state for the axion- 
modulus system (see Figure [T^ has a novel shape, with 
no oscillations, varying through behaviour like a cos- 
mological constant, radiation, matter, and quintessence 
as it begins rolling and goes through its various scal- 
ing stages. The total equation of state (Figure [T^ is 
again marginally perturbed by the presence of a signif- 
icant axion-modulus component near equality, and the 
final fractional density in the axion-modulus system is 
^ax+mod — 0.01. Compare to the previous example: the 
fractional density in the axion-modulus system is approx- 
imately the same (being fixed by the scaling solution by 
C), but makes a significant contribution at the slightly 
later time of a ~ 10~^, and correspondingly alters the 
total equation of state at this time. This demonstrates 
that there is control in the parameters over features in 
the total equation of state. The density at late times in 
this example is larger, since the tracking is never ended 
by axion oscialltions. 

Reducing the initial field for the axion in this example 
only changes the O's slightly, as one would expect when 
tracking is present, but moves the scale of the kink in Wtot 
to yet later times (not shown) . Increasing M undoes this 
change, as we expect since for slowly varying (p it is the 
X dependence of the potential that counts, and M'^cfy^ is 
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FIG. 16: Densities as a function of scale factor for the parame- 
ters of EguationfTilbut with the change M ->■ 10®, (pi 10~^. 
Notice that the axion field behaves completely as standard 
@, d]. The modulus begins a,t positive A like behaviour, before 
a fall into a negative potential dominated phase, signalling the 
onset of cosmological collapse. 




FIG. 17: Evolution of the modulus equation of state w as a 
function of scale factor for the parameters of Equation [141 but 
with the change M — >■ 10*, (f>i — > 10~^. The tracking solution 
is never quite found, as the equation of state makes a slow 
oscillation. When the modulus settles into its true negative 
potential minimum we see a the spike in the equation of state 
at a ~ 10"^ 



just a multiplicative factor in the exponential x potential, 
like B and D. It is the M'^cj)'^ term that is dominant for 
the parameters in question and as such the example is 
similar to the well studied case of a single scalar field 
with an exponential potential. 

Further variations now cause qualitative changes in the 
cosmology, which are demonstrated in Figures \W\ and [T71 
A higher mass M = 10* , and lower misalignment angle 
9i = 1 {4>i = 10~^) destroys the axion tracker and the 
axion behaves as standard, oscillating and making up 
a subleading fraction of dark matter, with the modulus 
potential dominated by the axion mass term and thus 
contributing a positive energy density. However, at later 
times the modulus falls into and bounces from its own 
potential minimum, eventually dominating with negative 
potential and again signalling rapid entry into a phase of 
cosmic contraction [52, [H^- During this evolution the 
modulus equation of state makes a slow oscillation as 
the field moves around in its potential and finds the true 
minimum. 

We finally note that contraction cannot be properly 
analysed in the framework presented here, since the 
Friedmann equation alone does not allow for it. The 
framework to use is the Friedmann acceleration equation 
for a Jss'], but applying that to this system, although 
simple in principle, is left for a future work. 



IV. CONCLUSIONS 

In this paper a model has been proposed that intro- 
duces a coupling between axions and moduli in the string 
axiverse. The coupling is motivated by the observation 
that the mechanism causing axion masses to distribute 
on a logarithmic scale, and thus for some to exist with 



masses in the range lO"'^"^ eV < < 10^^* eV where 
their cosmological dynamics can produce interesting phe- 
nomenology, is due to an exponential dependence of the 
mass on the size of cycles in the compact space. The 
sizes of these cycles are controlled by scalar fields called 
moduli, which are themselves dynamical. Moduli sta- 
bilisation is an important problem in string theory. We 
investigate the possibility that if there are light axions 
that roll in their potentials on timescales of cosmological 
interest then some moduli may also roll in their poten- 
tials, given by a general form for stabilisation, and that 
this vacuum destabilisation can be affected in turn by the 
presence of the axions. 

We then explore the consequences of having one cos- 
mologically relevant axion and allowing its counterpart 
modulus to also roll. The resulting system, in terms of 
the background expansion of the universe, is a simple 
one of two coupled scalar fields with a scalar potential 
containing a number of free parameters, decoupled from 
the other cosmic fluids. The potential causes the fields 
to have scaling solutions where the energy density tracks 
that of the dominant component. This destabilises the 
modulus, and the resulting evolution causes an evolution 
of the axion mass, altering the dynamics from the most 
simple case of a single decoupled axion. 

The axion-modulus system has a number of different 
phases in its evolution, of which we have identified some 
which may be of phenomenological interest. A common 
feature is that when tracking begins, the fraction of the 
total energy density in the axion-modulus system rises to 
an appreciable level and causes a non-standard evolution 
in the total equation of state in a manner similar to mod- 
els of Early Dark Energy. In the examples considered in 
this paper this resulted in a decrease of the equation of 
state by 0(10%) around the epoch of equality. 
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The fate of the axion-modulus system, and conse- 
quently the fate of cosmic expansion, then depends on 
the rate of decay of the axion mass and the amphtude of 
axion oscillations. If the axion oscillates and the ampli- 
tude of its oscillations decrease below some critical value 
then the modulus falls into its globally stabilising mini- 
mum at negative potential, which if this potential comes 
to dominate the energy density will signal the onset of an 
epoch of cosmic collapse. If the dynamical axion mass can 
be sufficiently damped by the modulus evolution then os- 
cillations cannot begin and the tracking solution remains 
stable. 

There is also the possibility that a modulus that looks 
to have been stabilised by Hubble friction at early times 
might be destabilised when its counterpart axion begins 
to roll at late times, and that this vacuum destabilisation 
may be observable indirectly through its effects on cou- 
plings of both gravity and the standard model, though 
the viability of this scenario in a realistic model requires 
much further attention. 

We have shown that these phenomena might be ex- 
pected as fairly generic since they are exhibited for values 
of the parameters spanning many orders of magnitude. 
We also note that in a scenario such as the string axi verse, 
where there is a plethora of light fields at the phenome- 
nologists disposal then it is possible to create cosmologies 
where phenomena like those described here may occur at 
multiple different epochs in the history and future of the 
universe, both by having multiple fields or by the two 
field dynamics spiralling towards some attractor and pe- 
riodically entering and exiting different regions in phase 
space. 

In the interests of simplicity no attempt was made to 
describe a universe where a significant and controllable 
fraction of the dark sector energy density is contributed 
by the axion-modulus system, except in the case where 
rapid decay of the modulus leads to cosmological collapse. 
This question of initial conditions will be the subject of 



a future work. 

Other future work will focus on delineating the regions 
of parameter space which give rise to the phenomena de- 
scribed in this paper and thus allow one to both construct 
desirable models for cosmological phenomenology based 
in string theory, and also to place limits on the possible 
values of parameters in any model of this type, which 
will be severely limited by the possibility of cosmological 
collapse. 

The coupling induces a tracking of energy densities 
in the axion-modulus system, which, just as it does for 
quintessence, may be fruitfully used to address problems 
of fine tuning, which are of a major concern for high fa 
axions. The tracking dynamics may also allow a modulus 
which appears naively to be stabilised in some negative 
energy density anti-de-Sitter minimum to in fact play 
the role of a cosmological constant leading to late time 
de-Sitter expansion, which may contribute in some way 
to the solution of the cosmological constant problem in 
string theory. 

This simple extension of the axiverse has a rich struc- 
ture and suggests models which we hope may be of use to 
both cosmologists and string theorists, and displays fea- 
tures which may, when fully investigated, be measured 
and constrained by cosmological experiments, particu- 
larly once the perturbations have been analysed to al- 
low computation of effects in the cosmic microwave back- 
ground and matter power spectrum. 
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